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In this paper we present a theoretical analysis of the performance of an 
m-phase coherent phase-shift keyed system in the presence of random 
gaussian noise and interference. An explicit expression is given for the 
probability of error of the phase angle of the received signal; we show that 
this probability of error can be expressed as a converging power series. 
We show that the coefficients of this series are expressible in terms of well- 
known and ivell-tabulated functions, and we give methods of evaluating 
the character error rates of the systems. We also show that this error rate 
is minimum when all the interference power is concentrated in a single 
interferer, and that it attains its maximum [P m ] mnx when the total inter- 
ference power is equally distributed amongst the K interferers. The limit- 
ing case when K goes to infinity is considered. The cases of K = 1, and 
m = 2, 4, 8, and 16 are treated in some detail, and the results are given 
graphically. The usefulness of the results presented in this paper is that 
the designer can have at his disposal very simple expressions with which to 
evaluate the performance of any given Coherent Phase-Shift Keyed sys- 
tem when the received signal is corrupted by both interference and random 
gaussian noise. 

I. INTRODUCTION 

The performance of coherent phase-shift keyed (CPSK) systems 
has been investigated by many authors; 1-5 in the transmission of in- 
formation the CPSK system has been shown to be one of the most 
efficient techniques for trading bandwidth for signal-to-noise ratio. 
However, the type of noise considered by these authors is almost 
always limited to be random gaussian noise although most authors 
admit that interference other than normal noise must be considered 
in the design of any modulation scheme for digital transmission. 

743 



744 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1969 

Consider the following situation. In the frequency bands above 10 
GHz where the signal attenuation resulting from rain storms could 
be very severe, close spacings of the repeaters are almost always 
mandatory for reliable communication from point to point and for 
all periods of time. 6 In such cases the problem of interference may be 
much more important than the problem of noise in the optimum detec- 
tion of the desired signal; hence it is very desirable to evaluate the per- 
formance of a CPSK system with co-channel and adjacent channel in- 
terference so that, for the selection of an optimum transmission scheme, 
comparative advantages of CPSK over other broadband modulation 
techniques (like FM) in combating interference can be determined. 

We consider in this paper the performance of a CPSK system when 
the received signal is corrupted by both interference and random 
gaussian noise.* We first discuss binary (2-phase) and quaternary 
(4-phase) CPSK systems and show that exact expressions can be 
obtained for their probability of error P m . These expressions are in 
the form of infinite power series which are shown to converge for all 
values of signal-to-noise ratio and for all signal-to-interference ratios 
above a certain level determined by the system. For m =. 2 and 4, 
these error rates are calculated and the results are given in graphical 
form. 

For m = 3 and for m > 4 we show that exact expressions for P m 
are very complicated functions of signal-to-noise ratio, and signal-to- 
interference ratios; in this paper we only indicate how these expressions 
can be obtained. However, we do obtain expressions for upper and lower 
bounds to P m and show that the difference between these two bounds is 
a monotonically decreasing function of signal-to-noise ratio, signal-to- 
interference ratios, and the number m of phases used in the system. 
For m ^ 4, signal-to-noise ratio p ^ 5 dB,f and for signal-to-interference 
ratio 1/L 2 ^ 20 dB, we show that this difference is less than 5 percent, 
and that the upper bound can be used as a good approximation to P„. 
For m = 8 and 16, we calculate these upper bounds and we present the 
results graphically. 

For a given amount of interference power, we show that the char- 
acter error rate is minimum when all the power is concentrated in 
a single interferer. If the total number of interferers is K we also 
show that the error rate P m reaches its maximum [P m ] m ax when the 
interference power is equally distributed among all the interferers. It 



*The word "noise" indicates random gaussian noise corrupting the desired 
received signal, 
t We use the notation b = a dB if 10 logioO = a. 
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follows that [P„,] max is a monotonically increasing function of K and 
attains its maximum when K goes to infinity. We show that the case 
of K going to infinity can be treated in a simple manner. 

For the computation of error rates P m (or upper bounds to P m , m 
> 2) it is necessary to calculate the central moments /i 2n 's of a certain 
random variable rj denned in terms of the K interfering carriers. For 
large values of K the conventional method of evaluating fj. 2n 's can be 
rather tedious; we give some simple methods of evaluating these mo- 
ments. 

In conclusion, this paper determines the performance of ?n-phase 
CPSK systems for the important case of signals corrupted by random 
gaussian noise and interference. The cases of m = 2, 4, 8, and 16 are 
treated in some detail. 

II. PHASE ANGLE DISTRIBUTION IN CPSK SYSTEMS 

Let us consider an ??i-phase CPSK system. We assume that there 
is a steady received signal* which is corrupted by random gaussian 
noise and interference. The gaussian noise is assumed to have zero 
mean and variance <r*. The signals under consideration consist of 
phase-modulation pulses of specified width transmitted at a known 
repetition rate ; we assume that there are K interferers, each interferer 
having the same form as the signal. 

If we assume that each signal transmitted has a duration T, the 
received signal waveform in the absence of noise during the iVth in- 
terval can be represented as 

s N (t) = (25)* cos (« * + 0), NT £ t £ (N + 1)T, (1) 

where S is the received signal power, w is the angular frequency of the 
signal, and will have some value in the discrete set 2irk/m, ^ k ^ 
m — 1, corresponding to the Nth message. All m messages are assumed 
to be equally likely. In the absence of noise and interference, the set of 
m possible received signals is described by a set of m equally-spaced 
vectors in the complex plane as shown in Fig. 1. The noise and inter- 
ference corrupting the signal distort the signal both in amplitude and in 
phase; a zero-phase signal (corresponding to A; = 0), as disturbed by 
noise and interference, is also shown in Fig. 1. 

If we now assume that powei in the jth interferer is /,-, the jth inter- 



* In this paper we do not consider the effects of fading on the error rates of 
CPSK systems. The effects of fading can usually be accounted for by a further 
integration of error rates obtained in this paper. 7 
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Fig. 1 — Phasor representation of CPSK signals for m = 4. 

ferer as received during the Nth interval can be represented as* 

i iN (t) = (27,)*cos{co,«-f-0, +»i\,NT£t£ (N + 1)T (2) 

where co y is the angular frequency of the jtb. interferer, 0/ is some value 
in the discrete set (2r/m) k, ^ k ^ m — 1, and the probability density 
Tr M (vi) of »i is given by 

0, otherwise. 

Since the K interferers are assumed to originate from K different 
sources, it is reasonable to assume that all /t/s are statistically in- 
dependent of each other and are also independent of gaussian noise 
n(t). 

The total received signal during the Nth interval can then be writ- 
ten as 

K 

r N (t) = (2S)* cos (« < + 0) + E (2/,)* cos («,/ + 0, + ft,) + n(0, 

NT ^ t g (JV + 1)T (4) 

where n{t) has zero mean and variance a 2 . 

Assuming that the receiver used in the system detects only the 
phase angle <£> of r N (t) and does not respond to its amplitude varia- 
tions,! we can write 8 



* We assume that all ijs's, I ^ j ^ K, are in the passband of the CPSK 
receiver used in the system. 

t This can be achieved in practice by using an ideal limiter at the front end 
of the receiver. If A(t)e ivW is the input to an ideal limiter, its output is given 
by A& iptn where A is a constant. 
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$ - tan" 1 ^ - u t (5) 

where t N (t) is the Hilbert transform of r N (t) and is given by 

m^-f^dr. (6) 

Let us write 

n(t) = I e cos (a t + 6) - I. sin (u t + 6). (7) 

We can show that I c and 7 a are two independent gaussian random 
variables each distributed with mean zero and variance <r 2 .* From 
(4) -(7) , we can now show that 

$ = d + tan" 1 

7. + E (2/,)* sin [(co, - co )i + 0,- - 6 + m,] 



j-i 



(2S)» + 7 C + E (27,)* cos [(«, - Uo )t +6,-6 + N ] 
i-i 

Let us now write 



(8) 



where 



and 



p = 


cr ' 




7. 

(25)* - 


y, 




I c 
(25)* " 


M, 




5 = 


E^ 

J = l 


sin X y , 


»7 = 


Efl, 

i-i 


cos X,- , 


«, 


-a 


I 



(9) 

(10) 
(11) 
(12) 
(13) 

(14) 

X,- = («, - Wo)« + 0, ~ 9 + M, . (15) 



* It is assumed that the spectrum of gaussian noise is symmetrical around the 
frequency u = wo . 
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Let us also denote the set {Xi , X 2 , ■ • • , X* , ■ • ■ > X*} of random 
variables X/s by X- 
We can now write eq. (8) as 

*> = 6 + tan" 1 = " + * (16) 

1 + U + 1} 

where 5 and i\ are functions of X. 

If K is a finite number, we can show 10 that the probability density 
■p v (rj) can be represented as* 

p,(i> = f P •"' n /•(«!) <«, ( i7 > 

where J (z) is the Bessel function of the first kind and of order zero. 
For K = 1, we can show that 10 

pm -■*&-&' M -*' as) 

0, otherwise. 

For K = 2, p v {rj) can be expressed in terms of elliptic functions, and 
for K > 2, no closed form expressions can be obtained for p„G?). In 
Ref. 10 Pr,(rj) has been expressed as a converging sum and has been 
evaluated for K = 10. It is easy to show that 

p,(„) =0, for | * | > S B, , ( 19 > 

j-i 

and 

f p,(i?)e _,,, di| = UMtRi). (20) 

J-00 1-1 

III. CPSK RECEIVER 

An ideal CPSK receiver is shown in Fig. 2. The ideal limiter re- 
moves all the amplitude variations of the received signal before it 
reaches the ideal phase detector of the system. We shall assume maxi- 
mum likelihood detection for our analysis of the receiver. Let us as- 
sume that the receiver shown in Fig. 2 has zero-width decision thres- 
holds as shown in Fig. 1. 



* We can also write similar expressions for p«(8) . 
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Fig. 2 — CPSK receiver. 

3.1 Error Rates for Binary CPSK Systems 

For a binary CPSK system the set of two possible received signals 
in the absence of noise and interference is shown in Fig. 3. The noise 
and interference corrupting the desired signal distort the signal both 
in amplitude and in phase; a zero-phase signal (corresponding to 
k = 0) as disturbed by noise and interference is also shown in Fig. 3. 

When the message k = is sent, and when the phase angle $ of the 
received signal lies in the second and third quadrants of the complex 
plane shown in Fig. 3, an error is made in detecting the received signal. 
For a given p 2 , and for an arbitrary set of X,'s let us assume that the 
origin of the gaussian noise vector is at point G in Fig. 3. When the 
terminus or tip of the gaussian noise vector lies in the left half of 
the complex plane (the shaded portion of Fig. 3) an error is made by the 
receiver. Since I e and I. are two independent gaussian random vari- 
ables and since they are distributed independently of X/s, the prob- 
ability P 2 (X) that the terminus of the gaussian noise vector lies in the 




Fig. 3 — Phasor representation of CPSK signals for m = 2. /« and I. are the 
in-phase and quadrature components of gaussian noise corrupting the desired 
received signal. 
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left half of the complex plane is given by* 
P 2 (l) = Pr[-oo < I, < oo] 

•Pr [-00 < I e < -{{2S) h + E (2/,)* cos X,)j 

■j (.-) (25) »+ 2Z (2/-/)»ooa\/ I 

"flk/J '"' /»P(-'/V)*. (2D 

We can show from Equation (21) that 

P 2 (X) = | erfc [p + p*], (22) 



where 



erf (x) = 4 f exp (-w 2 ) dw (23) 



and 

erfc (x) = 1 - erf (x). (24) 

The character error rate P 2 for a binary CPSK system is, therefore, 
given by 

P 2 = E[P 2 (\)], (25) 

where E[P 2 (\)] represents the mathematical expectation of the random 
function P 2 (h)- 

From Equations (22) and (25) we have 

P 2 = i£[erfc {p + pr?}]. (26) 

We now note that we can write 11 ' 12 

erfc [x + «] = erfc [x] + h exp (-z 2 ) £ (-l)*H,_,(af) jt , (27) 

where #„ (.t) represents the Hermite polynomial of order n. The series 
converges for all values of x + z such that 

x + z ^ 0. (28) 

From Equations (26) and (27) we have 

P* = I erfc (p) + -\ exp (-p 2 ) £ (-l)^^(p) £ S(0- (29) 



*The notation Pr[a < a; < o] denotes the probability that the random vari- 
able x satisfies the inequality a < x < b. It may also be noted that P 2 (X) is a 
conditional probability conditioned on X. 
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Let us denote by ^ the nth central moment of 17. 8 It can then be 
shown that 10 

Ihi+i = 0, 1 = 0,1,2,"-. (30) 

We, therefore, have 

P 2 = I erfc (p) + ^ exp (-p 2 ) g tf^Cp) ^j p,, . (31) 

The series given in Equation (31) converges for all values of p and 
Rfs such that 

P + pv ^ for all &. (32) 

From Equations (13) and (32) we can show that the series converges 
when 

0^1, (33) 

where 

Q = T,Ri- (34) 

j-i 

Equation (34) states that the sum of the normalized amplitudes of all 

the interfering carriers may not exceed the normalized amplitude of 

the desired signal. This is not a very stringent requirement and it is 

almost always satisfied when low error rates are desired. 

Since we also know that* 

j £«,}* gar, (35) 

when Equation (33) is satisfied, we have 

The expression S/Ij denotes the signal-to-interference ratio of the 
;'th interfering carrier. 

When there is only one interfering carrier we can show that, 



(2Q ! 



*« - RU S . (37) 



♦Equation (35) states that the arithmetic mean of a set of real variables 
is always greater than or equal to its geometric mean. 
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and equation (31) can be written as 



[<r 



P 2 = h erfc ( P ) + ~j exp (-p 2 ) g ffM-i(p) **ggr- ( 38 ) 

The series in equation (38) converges for all signal-to-interference 
ratios such that 

1/8 ^ 1. (39) 

The values of P 2 have been calculated from equation (38) and the 
results are given in graphical form in Fig. 4.* 

Notice that we need to calculate only the even order moments 
fi 2 »'s of the random variable v in determining P 2 from equation (31). 
Some methods of calculating these moments are given in Appendix A. 

3.2 Error Rates for Quaternary CPSK Systems 

Let us now consider a 4-phase CPSK system. For this system the 
set of four possible signal phasors and the four optimum decision 
thresholds are shown in Fig. 5. A signal phasor (corresponding to 
k = 1) as disturbed by noise and interference is also shown in Fig. 5. 

For a given set of A/s let us assume that the gaussian noise is rep- 
resented by a vector from the point G. If the message k = 1 is trans- 
mitted, an error is made if the received phase angle lies in areas 
marked 1, 2, and 3. The phase angle of the received signal will lie 
in areas marked 1, 2, and 3 if the terminus of the gaussian noise 
vector lies in this area of the plane. 14 

We notice that 

GA = (2S) i sin | + £j (21,)* sin (| + X,) , (40) 

and 

GB = (25)* sin | + fj (27,)* cos (| + h) ■ (41) 

Let us denote by n* llM , ■ • ■ , *„(A) the probability that the terminus 
of the gaussian noise vector lies in area 

Om 

* The results obtained in Fig. 4 indicate that the error rates obtained in Refs. 
13, 14, and 15 agree well with those obtained in this paper. 

tThe notation U"_iA\- denotes the union of all elements of the set{«i,«2,* ■ •,*■}. 
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Fig. 4 — Error rates for a 2-phase CPSK system with one interferer. 
We can show from Fig. 5 that 



n ll2 (X) = i erfc I psin| + p S #, cos (| + X,jJ , 
n 2 . 3 (X) = \ erfc I p sin | + p £ R, sin (| + A, J J , 



(42) 
(43) 



and 



n 2 (X) = \ erfc 



|_psin| 



- + p 2 Ri cos 



5 + *)] 



• erfc [p sin | + p £ fi, sin (| + X,) J • (44) 

The probability P 4 (X) of an error due to noise is, therefore, given by 

P&) = n,. 2 (x) + n 2 .,(x) - n 2 (x). (45) 

The probability of an error due to noise and interference is therefore 
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QUADRATURE AXIS 




Fig. 5 — Phasor representation of CPSK signals for m = 4. 7 U and 7, are two 
orthogonal components of gaussian noise. 

given by 

P 4 = E[PMl (46) 

From equations (27), and (42) through (46) we can show that 



P 4 = erfc 



\ P sin | J - i erfc 2 p sin | J 



+ -rsi exp y- P 2 sin 2 |JJ2 - erfc |^p sin |J| 
' 2 H«-i(p sin |) J^ *, - ; exp (-2p 2 sin 2 1) 






P MaVa* 



(2£)!(2;)! 



(47) 



where m*/ i2 ,'s are given by 



'.« = 7T^ de i dd * * ' ' i d M 23 «/ cos 0,- , 

(2t) ■'o ■'o •'o ki— l J 

.</£B t 8ia0tj • (48) 
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For a given set of JR/s, m*< i2 ,'s may be evaluated from equation (48). 
For K = 1, we can show that 

(20! (2s)! 



2U + .) 



..* _ D2(<+« 

M2<.2i — I* 



2 2{l+,) t\s\{l + s)\ 



(49) 



For K = 1, we have calculated P 4 from equation (47) and the 
results are presented in Fig. 6. 

We can again show that the series given in equation (47) con- 
verges for all values of p and R/s such that 



o ^ ■ * l 



(50) 



For K = 1, equation (50) becomes 

S/I ^ 2. (51) 

Equation (50) is usually satisfied by systems encountered in practice. 




10 II 
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Fig. 6 — Error rates for a 4-phasc CPSK system with one interferer. 
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3.3 Error Rates ]or Multilevel CPSK Systems 

In this section we shall investigate the performance of a multilevel 
(m ^ 3) CPSK Systems and indicate a method in which an exact 
expression can be obtained for the probability of error of the system. 
This exact expression is a very complicated function of signal-to-noise 
ratio and fl/s; we do not obtain this expression in this paper. However, 
we obtain upper and lower bounds to P m and show that the difference 
between these two bounds is a monotonically decreasing function of 
p, m, and signal-to-interference ratios. For K = 1, m ^ 4, p 2 ^ 5 dB, and 
8/1 ^ 20 dB, we show that this difference is less than 5 percent of 
the lower bound, and hence the upper bound is a good approximation 
to P m when low error rates are desired. 

A signal phasor corresponding to k = as disturbed by noise and 
interference is shown in Fig. 7. For a given set of X,-'s let us again assume 
that random gaussian noise is represented by a vector from the point 
G shown in Fig. 7. If the message A; = is transmitted, an error is 
made if the terminus of the noise vector lies in areas marked 1, 2, and 3. 

We can show that * 



n,. 2 (X) = i erfc J^psin ^ + p E R, sin \± - x,jj 



(52) 



and 



n 2 ,(X) = \ erfc [ P sin i + p ± R, sin (^ + X,)] • (53) 

The probability of error due to noise is, therefore, given by 

PM = n ll2 (x) + n 2i3 (x) - n 2 (x). (54) 

By looking at Fig. 7 we can see that no simple expression can be 
obtained for n 2 (X) (except when m = 4). U 2 (W denotes the probability 
that the terminus of the gaussian noise vector lies in area 2; we shall 
now obtain upper and lower bounds to n 2 (X). Assume that 

psin-- ptiJ^O (55) 



Note that 



and 



GA = (2S)" sin i + t (2//)' sin (^ - X,-) 
GB = (2S)> sin ^ + t (2I*) 1 sin (^ + X,-)- 
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Fig. 7 — Phasor representation of CPSK signals for m = 8. 7 U and h are the 
in-phase and quadrature components of gaussian noise. 

so that n 1>2 (X) and n 2i3 (X) are nonnegative for all values of X. If equa- 
tion (55) is satisfied, it is easy to see (see Fig. 7) that 

n 9 (fc) ^ for all X, (56) 

and n 2 (X) reaches its maximum when* 



n = -T,R,- = -a. 



(57) 



For this value of rj it can be shown (see Fig. 8) that 
n 2 (X) =-K I "" exp (-y 2 /2a 2 ) dy f exp (-x 2 /2a 2 ) dx 

7TO" J -„ Jo 



or 

n 2 (X) - -i> f " exp (-y 2 /2a 2 ) dy f " exp (-x 2 /2a 2 ) dx (58) 

ir<r J u , Jo 

where 

y = (2S)*[1 - Q]. (59) 

Since we always have 

^ exp (-x 2 /2ff 2 ) ^ 1 for all real x, (60) 



* We can show that j? = — fi when all X/s are odd multiples of jr. 
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't 




Fig. 8 — Computation of lower bound to P m . 
we have 

n 2 (X) g t -^4^ f (y - Vo) exp (-y 2 /2a>) dy. 

Equation (61) can be simplified to* 

„,,, ^ r. tan7r/w r 2,-, n \2-, 
H 2 (X) ^ Q m0 = — exp [-p (1 - Q) J 

IT 

•[1 - (7r)V(l - 0) exp [p 2 (l - Q) 2 ] erfc {p(l - 0)}]. 
From equations (54) , (56) , and (62) we have 

n ll2 (x) + n li3 (x) - Q m0 ^ PM £ n 1>a (x) + n lf8 (x). 

Since 

P m = J[P.a)L 
we can show from equations (52), (53), (55), and (63) that 

Q m - Q m0 ^ P m ^Q m 
where 

Q m = erfc (p sin ^ 

, 2 / 2 ■ 8T \A g 4 8k m) 2< 

+ ^exp^-p sm -J g ^ p *,. 



(61) 



(62) 
(63) 
(64) 
(65) 



(66) 



For large values of p and small values of fi, Q m o is approximately equal to 

tan 7r/w exp [— p 2 (l — Q) 2 ] 
2tt p 2 (l - G) 2 
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The series given in equation (66) converges if equation (55) is 
satisfied, or if 

fi<sin-- (67) 

m 

When low error rates are desired, equation (67) must be satisfied. 
Equation (65) gives an upper and a lower bound to P m ; as can be 
seen from equation (62) the difference Q m0 between these two bounds 
is a monotonically decreasing function of p, m, and signal-to-inter- 
ference ratios. From equation (65) we have 

< Pm ~ Qm ^ 0. (68) 



Q m - Q m0 = P 

For K = 1, #, = tV> and for m = 4, 8, and 16, we have plotted in 
Fig. 9 Q m0 /(Q m - Qmo) as a function of p 2 . From Fig. 9 we see that 
Qmo/(Qm — Qmo) is less than 5 percent for p 2 ^ 5 dB and for m ^ 4. 
We can, therefore, use Q m as a good approximation to P m for high values 
of signal-to-noise ratio (p 2 ^ 5 dB) and for high values of signal-to- 
interference ratio (l/Ri ^ 10 dB). 

In these cases we then have 

P m tt erf c ( p sin — J 
\ ml 

2 / \ " H2k ' 1 \ psm m) 
+ ^ exp (-p 2 sin 2 f n ) g m] p 2 W • (69) 

For K = 1, and for m = 8 and 16, the values of P m obtained from 
equation (69) are given in Figs. 10 and 11. The error made in this 
approximation can be estimated from equation (68) . 

IV. ERROR RATE AS A FUNCTION OF NUMBER OF INTERFERERS 

Let us now investigate how P m varies as a function of K for a total 
given interference power. Let us assume that the total interference 
power is some number SL 2 where 

E /,• = SL 2 . (70) 

This power SL 2 can be distributed among the K interferers in a 
variety of ways; every one of these distributions will in general lead 
to a different character error rate of the system. Let us find out those 
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SIGNAL-TO-NOISE RATIO IN DECIBELS 

Fig. 9 — Qmo/(Qm — Qmo) as a function of p. 

distributions of power (if they exist) which make this character error 
rate a maximum or a minimum. 

4.1 Error Rates for K Inter jerers 

Let us first consider the case when p ~2> 1 and ft « 1. In this case the 
series corresponding to P m (or Q m ) converges very rapidly; let us say 
that the first N terms of the series are sufficient to evaluate P m to the 
desired degree of accuracy. 

For all I and z, we have 



#2,-1(2) = 2zH 2l - 2 (z) -2{2C- 2)# 2 ,_ 3 (2). 
From equation (71) we can show that 



H 2i ^[psm^j ^ 0, 1 gjgN, 



(71) 



(72) 
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if 



P ^ 



2N - f 

. 7T 

sin — 
m 



.V > 1. 



(73) 



If Equations (72) and (73) are satisfied notice from Equations (31) 
and (69) that P m 's are monotonically increasing functions of n 2 i's, 
I 2> 1. For a given n 2 , it can be shown from Equation (13) that ^ 2 ,'s 
I Si 2, reach their minimum when ft is minimum and they reach their 
maximum when ft is maximum. 

We can then say that P m 's (or Q ra 's) attain their minimum when ft 
is minimum and that they are at their maximum when ft is maximum. 

From Figs. 3, 5, and 7 this seems to be true for all values of p and ft 
which satisfy Equation (55). 

Let us now find out when ft is minimum for a given value of signal- 
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Fig. 10 — Error rates for an 8-phase CPSK system with one interferes 
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Fig. 11 — Error rates for a 16-phase CPSK system with one interferer. 
to-interference ratio. The signal-to-interference ratio 1/L is given by 



K r 

Clearly fi is minimum when 

/, = SL 2 , l^j^K, 
and 



(74) 



(75) 



(76) 



We can then say that the character error rate P m is minimum when the 
total interference power is concentrated in a single interferer. 
Now from equations (14), (34) and (74), O is a maximum when 



J) 
dl 



;[£(?-<£©] = <>• **'**- 



(77) 
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e is a constant and is the Lagrange multiplier used in finding the 
extremum of o. 

Solving equation (77) we observe that O is a maximum or that 
P„, is a maximum when 

I, = ^f , l^j^K (78) 

or that the total interference power is equally distributed among the 
K interferers. 

Let us now assume that K is a variable number. It is clear from 
equation (78) that [P,„] ma x is a monotonically increasing function of 
K. 

4.2 Error Rates for a Large Number of Interferers* 

Let us now consider the limiting case when K goes to infinity and 

f:i t -8L\ (79) 

We can show 10 that the probability distribution function off 

y(t) = £j (21,)* cos {«,* + 0, + p,} (80) 

j-i 

as K goes to infinity approaches that of gaussian noise with mean 

zero and variance SL 2 under certain conditions. 

In this case we have from equation (4) 

r N {t) = (2S) 1 cos ( Uo < + d) + y(t) + n{t). (81) 

Since y(t) and n{t) are independent gaussian random variables their 
sum 

b(t) = y(t) + n(t) (82) 

is also a random gaussian variable with mean zero and variance SL 2 

+ <T 2 . 

From equations (81) and (82) we can write 

r„(t) = (25)* cos (o> * + d) + b(t) (83) 

where b (t) is a gaussian random variable. 



* The results of this section are applicable for any signal-to-noise ratio and 
any signal-to-interfcrence ratio. 

t Ruthroff has shown that for K & 50 the distribution of y(t) can be considered 
(o lie gaussian in practice for the computation of distortion in PM systems. 10 



764 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1969 

The case where r N (£) can be described by equation (83) has been 
considered in detail in Ref. 17;* we can easily determine the deteriora- 
tion in performance produced by interference from the results pre- 
sented in that paper. For example, suppose that m = 4, S/a 3 = 16 dB, 
and U = -16 dB. Clearly 

£ + L 2 = -13dB (84) 

and P 4 from Ref. 17 is given by 

P 4 = 7.9 X 10" 6 . (85) 

For the calculation of the effect of interference in CPSK systems, 
we note that we have not shown the validity of the gaussian approxi- 
mation of i/(0 for K » 1. However, this assumption seems to be 
justified for large signal-to-noise ratios and small interference-to- 
signal ratios. 15 

In conclusion, this section gives methods of evaluating character 
error rates of CPSK systems for all values of m and for all values of 
K. It shows that the error rate P m is minimum when all the inter- 
ference power is concentrated in a single interferer and that it attains 
its maximum value [P„,] m ax when the interference power is equally 
distributed amongst all the interferes. We further show that [P„,] mn x 
is a monotonically increasing function of the number K of interferers. 
We also show that the case, K going to infinity, can be treated and that 
the deterioration in performance produced by interference can be de- 
termined. 

V. CONCLUSIONS 

A method to evaluate the character error rates of CPSK systems 
has been presented in this paper. The received signal is assumed to be 
corrupted by both interference and random gaussian noise. When the 
number of interferers is very large it can be shown that the inter- 
ference and random gaussian noise can be combined together to give 
rise to an equivalent noise source having gaussian properties. The 
variance of this random variable is the sum of variance of random 
gaussian noise and total interference power. In this case the analysis 
of the CPSK system can be done by methods presented in Ref. 17. 

When K is a finite number and when m = 2 or 4, exact expressions 



* The results presented in this paper for S/I = oo are also sufficient to deter- 
mine P m for a large number of interferers. 
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are given for the probability of error P m . When m ^ 3, upper and lower 
bounds to P m are derived. We show that the difference between these 
two bounds is a monotonically decreasing function of signal-to-noise 
ratio p 2 , signal-to-interference ratio 1/L 2 , and the number m of phases 
used in the system. For K = 1, m ^ 4, p 2 ^ 5 dB, and l/R, ^ 10 dB 
we show that this difference is less than 5 per cent, and that the upper 
bound can be used as a good approximation to P m . 

We then show that for any ra-phase CPSK system the character 
error rates can be expressed in terms of the central moments of a 
certain random variable -q and that they can be calculated to any 
desired degree of accuracy by using a set of tables or by using a 
digital computer. 

For a total given interference power we show that the character 
error rate P m attains its minimum when all the power is concentrated 
in a single interferer, and that it reaches its maximum [P ni ] ma x when 
the power is equally distributed among all the K interferers. It is also 
shown that [P m ] nmx is a monotonically increasing function of K. 

The cases of K = 1, m = 2, 4, 8, and 16, have been treated in some 
detail and the results are given in graphical form. The required signal- 
to-noise ratio for any value of signal-to-interference ratio can be 
determined from these figures. 

The usefulness of the presented results is that they provide the 
designer with some relatively simple expressions with which to eval- 
uate the performance of any given CPSK system with interference 
and random gaussian noise. The only quantities he must have at his 
disposal are the central moments of a certain random variable rj 
denned in terms of the K interfering carriers. 
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APPENDIX 

Evaluation of Central Moments of rj 

In the computation of character error rates for CPSK systems it is 
necessary to calculate the even order moments of the random variable 
77; we shall give in this section two alternate methods to evaluate these 
moments. 
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By definition pan is given by 
to* = 7T7£ / de i \ dd 2 "- de K \ £ R, cos 6, . 

\2ir) J o Jo J o Lj = i J 

By the multinomial theorem 

[k "]2n (2n)' K 

z r, cos oA = z k- n aw cos-' *, 

y-i 
where n/s are positive integers such that 

K 

2 W I = 2W - 



(86) 



(87) 



(88) 



Since 0/s are statistically independent of each other, and since 
M2<+i = for all I, we have from Equations (86) and (87)* 



^-E-^nw- (w ' )! 



Un ( l "* 



2 «< 



rtiu 



(89) 



where n/s are a set of even positive integers satisfying Equation (88). 

Even though equation (89) gives an exact expression to evaluate 
fi2n's, it can be rather tedious to evaluate fi 2n 's from equation (89) 
for large values of n and K. We shall therefore give an alternate 
method to evaluate the central moments of the random variable -q. 

It can be shown that the probability density function p v (v) of the 
random variable rj can be expressed as 10 



*M - 2fi 



11 = 1 " 1 = 1 






The 2nth moment of 77 can be represented as 

M2n = / 3 a "p,(z) dz. 

J-a 

From equations (90) and (91) we can show that 



(90) 



(91) 



n^ 






E(-D* 



(2n)! 



[2n - 2k + 1] ! (£r) s 



(92) 



* For X = 1, equation (89) reduces to equation (37). 
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It can be seen that the infinite series appearing in equation (92) 
converges rapidly for all values of R/s; we need take only a finite 
number of terms from equation (92) to estimate p^'s. It is, therefore, 
easier to evaluate /A2»'s from equation (92) than from equation (89) 
when there are a large number of interferers, and we have to take a 
large number of terms in estimating P m . 
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